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I derive formulas for the electrostatic potential of a charge in or near a membrane modeled as one 
or more dielectric slabs lying between two semi-infinite dielectrics. One can use these formulas in 
Monte Carlo codes to compute the distribution of ions near cell membranes more accurately than 
by using Poisson-Boltzmann theory or its linearized version. Here I use them to discuss the electric 
field of a uniformly charged membrane, the image charges of an ion, the distribution of salt ions 
near a charged membrane, the energy of a zwitterion near a lipid slab, and the effect of including 
the phosphate head groups as thin layers of high electric permittivity. 



I. CELL MEMBRANES 

The plasma membrane of an animal cell and the mem- 
branes of the endoplasmic reticulum, the Golgi appara- 
tus, the endosomes, and other membrane-enclosed or- 
ganelles are lipid bilayers about 5-nm thick studded with 
proteins. The lipid constituents arc mainly phospho- 
lipids, sterols, and glycolipids. 

Of the four main phospholipids in membranes, 
three — phosphatidylcthanolamine (PE), phosphatidyl- 
choline (PC), and sphingomyelin (SM) — are neutral, and 
one, phosphatidylserine (PS), is negatively charged. In a 
living cell, PE and PS are mostly in the cytosolic layer 
of the plasma membrane; PC and SM are mostly in the 
outer layer [ij; and the electrostatic potential of the cy- 
tosol is 20 to 120 mV lower than that of the extracellular 
environment. 

After pioneering work by Gouy, Chapman, and Wag- 
ner 0, and by Onsager and Samaras 0], many scien- 
tists have studied the electrical properties of cell mem- 
branes 0, • This paper presents the exact electrostatic 
potential due to a charge in or near a membrane in the 
continuum limit in which the membrane is taken to be 
one or more dielectric slabs lying between two different 
infinite dielectric media. Because of the superposition 
principle, this monopole potential also gives the multi- 
pole potential due to any array of charges in or near a 
membrane. 

One can use these formulas to simulate the interactions 
of ions with other ions and with fixed charges near mem- 
branes while modeling water and other neutral molecules 
as bulk media. For instance, one can use them in Monte 
Carlo simulations to compute the behavior of salt ions 
and protons in water near neutral or charged membranes 
even in the presence of fixed charges of arbitrary geom- 
etry. This method is more accurate than the Poisson- 
Boltzmann mean-field approximation and much more ac- 
curate than its linearized version Q . These formulas also 
provide a context for and a check on all-atom computer 



simulations [1, 0] . 

As early as 1924, Wagner [2| noted that an ion in wa- 
ter near a lipid slab induces image charges that repel 
the ion. No mean-field theory can describe this simple 
effect. But work-arounds are available for the Poisson- 
Boltzmann theory Q. 

Formulas for the electrostatic potential of a charge in or 
near a cell membrane modeled as a single slab are derived 
in section ini As pedagogical illustrations of their utility, 
I use them to compute the electric field of a charged mem- 
brane in section IIIII and the response of bound charge to 
an ion in section IIVI In section |Vl I use them to sim- 
ulate the distribution of salt ions near a charged mem- 
brane. I discuss the Debye layer in section IVII and the 
energy of a zwitterion near a lipid slab in section I VIII 
In section rVlIIl I calculate the potential of a charge near 
a membrane modeled as several dielectric layers of dif- 
ferent permittivities between two different semi-infinite 
dielectrics. I use this analysis in section IIXI to model 
a phospholipid bilayer as a lipid layer bounded by two 
layers of head groups of high electric permittivity. The 
phosphate head groups cause the membrane to attract 
rather than to repel ions. I summarize the paper in sec- 
tion H 



II. THE POTENTIAL OF A CHARGE IN OR 
NEAR A LIPID BILAYER 

In electrostatic problems, Maxwell's equations reduce 
to Gauss's law V • D = which relates the divergence 
of the electric displacement D to the free-charge density 
Pf (not including the polarization of the medium), and 
the static form of Faraday's law V X E = which implies 
that the electric field E is the gradient of an electrostatic 
potential E = -~W. 

Across an interface with normal vector fi between two 
dielectrics, the tangential component of the electric field 
is continuous 



n X (E2 - El) = 



(1) 
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while the normal component of the electric displacement 
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jumps by the surface density a of free charge 



fi • (D2 - Di) = cr. 



(2) 



In a hnear dielectric, the electric displacement D is the 
electric field scaled by the permittivity e of the material 
D = eE. 

The lipid bilayer is taken to be flat, extending to in- 
finity in the x-y plane, and of a thickness i w 5 nm. The 
interface between the extracellular salty water and the 
lipid bilayer is at 2; = 0. The permittivity of the lipid 
bilayer is about twice that of the vacuum k, 2eo; those 
of the extracellular environment e.^, and of the cytosol ec 
are about 80 times eo- 

The potential of a charge g at a point (0,0, h) on the 
z-axis is cylindrically symmetric, and so Bessel func- 
tions are useful here. In cylindrical coordinates with 
p = a/x^ + y^, the functions J„(fcp) e™*^ e^*^^ form a 
complete set of solutions of Laplace's equation, but due 
to the azimuthal symmetry, we only need the n = func- 
tions Jo{kp) e^^^. We will use them and the relation Q 
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^/p^ + {z- hy Jq 



dk Jq {kp) e 



— k\z — h\ 



(3) 



to represent the potential of point charge at (0,0, h). 

If the charge q is at (0,0, /i) in the water above the 
membrane (/i > 0), then we may write the potentials in 
the extracellular water V.^ , in the lipid membrane , 
and in the cytosol as 



g ^-k\z-h\ 



47re,, 



FJ(p,z)= / dkUkp) 
Jo 

V^{p,z)^ / dkJo[kp)[m{k)e^'^f[k)e 
Jo 

poo 

kz 



u{k) e 



— kz 



— kz 



y-(p,z) 



dk Joikp) d{k) 



(4) 







Imposing the constraints (IT]&[2]), writing u(k) as u, mik) 
as m, and so forth, and setting /? = qe''^'^ /AnCui and 
y = e^*^*, we get the four equations 

m + f — u ~ f3 
egm - egf + e,„u = e„,/3 
eim - eiyf - = 
m + yf — d = 0. 

In terms of the abbreviations 
1 



(5) 



P 



(e™ + e^) and = - (cc + e^) 



and p' 
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Cc + ee 



(6) 



their solutions are 



u{k) ^ p 



p - p'/y 

l-pp'/y 

m{k) = /3 



d(k) = /3 



suje 1 - pp'/y 
eto p'/y 



f-wi 1 - pp'/y 



f-w(.f-c(. 1 - pp'/y 



(7) 



Inserting these solutions into the Bessel expansions ^ 
for the potentials, expanding their denominators 



1 - pp'/y 



-2nkt 



(8) 



and using the integral ([3]), we find that the potential 
V^'{p,z) in the extracellular water due to a charge q at 
(0, 0, h) in that water is 



Kr(p,^) 



47re„ 



/ + iz + hf 

{pp'r-^ 



(9) 



p'(i-/)E 



1 + (z + 277.f + hy 



in which r = -^/p^ + (z — /i)^ is the distance to the charge 
5, and p'^ + (z + h)^ is the distance to the principal 
image charge pq. Similarly, the potential V^™ in the lipid 
bilayer is 



vrip,z) 



_^ _ 2nt - hf 



^p-^ + {z + 2{n + l)t + hY ^ 
and the potential V^^ in the cytosol is 

qu ^ {pp'T 



(10) 



V:{.P.z)^ 



47re, 



£erf,^g + (z - 2ni - 



(11) 



If the charge 5 is in the lipid bilayer at (0, 0, h) with 
— i < ft- < 0, then the Bessel representations of the po- 
tentials are 



vi{p,z) 

V,\p.z) 



dk Ja{kp) u{k) e 



-kz 



dk Ja{kp) 



_!i_-k\z^h\ 



+m{k) e^'- + f{k) e-'^^ 

D 

dkJo{kp)d{k) e'=^ 



(12) 
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With 7 = qe''''' / 4:TTe£ and x = e '^^^ , the Maxwell con- 
straints ([1] & [5]) give us the equations 

u — m — / = 7 
e^u + e^m - ee f = e^j 
d^m. — yf = xj 
ecd - et/m - eiy f ^ eixj. (13) 

whose solutions are 



u{k) 
m{k) 

f{k) 
d{k) 



f-tP 1 - p'x/y 
e^,t 1 - pp' Iv 

1 - pp'/y 
p' 13 X — p 



y I -pp'/y 
eiP X - p 



(14) 



£ct 1 - pp'/y' 

After putting these solutions into the Bessel expansions 
([T2|) and using the integral (|3]) and the denominator sum 
([5]), one finds that the potential in the extracellular 
water due to a charge q at (0, 0, h) in the lipid bilayer is 



E 



ipp') 



E- 



^^p2 + (z + 2nt-h)^ 

p' [ppT 



The potential V/ in the lipid bilayer is 



(15) 



y/(p,z) 



47ref 



[pp') 



l\\n\ 



p' {pp'r 



(16) 
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and that in the cytosol is 



V^ip^z) 



47rec 



(pp')" 



y— 

+ - 2nt - /i)2 



p (ppQ" 

+ (z - 2nt + /i)2 



(17) 



Finally and somewhat redundantly, we turn to the case 
of a charge q in the cytosol at (0, 0, h) with h < —t. Now 
the Bessel expansions of the potentials are 



V^{p,z)= dkJo{kp)u{k)e' 
Jo 



-kz 



(18) 



V[{p, z)= dk Jo{kp) [m{k) e^' + /(fc) e 
Jo 

F/(p,z)= / dkJoikp) 



d{k)e'' 



The continuity conditions ([T] & [2]) give us the four equa- 
tions 



u — m — f — 
e^u + e^m - e^/ = 

d - m -yf = - /3y 
Ccd - etm + eiyf = I3cy 



(19) 



whose solutions are 
u{k) 

m{k) 

fik) 

d{k) 



f-wif-ct 1- pp'/y 
_ ^ pP 
ect 1 - pp'/y 

eel 1 - pp'/y 



py-p 

1 - pp'/y ■ 



(20) 



Thus using ([3] & [8]) in , we find that the potential 
in the extracellular water due to a charge q at (0, 0, h) 
in the cytosol is 



v:{p.z) 



■E 



{pp'Y 



47r ey,i eel ^ + (z + 2nt - h)"^ 
The potential V[ in the lipid bilayer is 



(21) 



V[{p,z) 



q 



Ait ece 



E 



{pp') 



'E' 



^ + (2; - /i + 27lt) 

(pp')« 



+ + /j _ 2„t)2 

The potential in the cytosol is 



(22) 



K^(p,z) 



1 



Aire 

p{i-p'')y: 



+ {z- 2nt + hY 



(23) 



n=0 



in which r is the distance to the charge and 
y/fP + (z + h + 2FY is the distance to the principal im- 
age charge p'q. 

Inasmuch as 1—p'^ — ^w^e/^wf series ((MTT|) for the 
potentials V^{p,z), V^{p,z), and V^{p,z) agree with 
those derived by the method of image charges 0. The 
first eight terms of the infinite series (|9l-[TTl [TSl— [T7l & 
[2TI— [23]) give the potentials to within a percent. They 
are fast. 

The first 1000 terms of the series P [TUl fc [TT|) for 
the potentials V^'"(p,z), V//'{p,z), and V^{p,z) (right 
curve, red) and (I1T1[22 &[231) for the potentials Vf{p, z), 
V^{p,z), and V/^{p,z) (left curve, blue) are plotted in 
Fig. [T] (V) for p = 1 nm as a function of the height z 
(nm) above the phospholipid bilayer for a unit charge 
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FIG. 1. (Color online) The electric potentials V^(p, z), 
Vr(P,^), and V^{p,z) dllTOlinil and VZ{p,z), Vl{p,z), and 
V^{p, z) (l2T][22l[2ll) are plotted (V) for p = 1 as a function of 
the height z above the phospholipid bilayer for a unit charge 
q — \e\ in the cytosol at (p, z) — (0, —6) (left curve, blue) and 
in the extracellular salty water at (0, 1) (right curve, red). 
The lipid bilayer is between z = —5 and z = 0; distances are 
in nm; and the permittivities are ew = = 80eo and ei — 2eo. 
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FIG. 2. (Color online) The electric potentials V^{p,z), 
V/(p,2), and Vi(p,2) (UllIlSl (V) for p = 1 as a func- 

tion of the height z above the phospholipid bilayer for a unit 
charge q = |e| in the phospholipid bilayer at (p, z) = (0, —4.5) 
(left curve, blue), (0,-2.5) (middle curve, magenta), and 
(0,-0.5) (right curve, red). Distances and permittivities are 
as in Fig. [T] 



q = \e\ in the extracellular medium at (p, z) = (0, 1) nm 
(right curve, red) and in the cytosol at (p, z) = (0, — 6) 
nm (left curve, blue). The errors due to the truncation 
of these series at 1000 terms are less than 1 part in 10^^. 

The first 1000 terms of the scries dm [HI & [HI) for the 
potentials Vf{p,z), V^{p,z), and V^{p,z) are plotted in 
Fig. [5] (V) for p = 1 nm as a function of the height z 
(nm) above the phospholipid bilayer for a unit charge 
q ~ \e\ in the bilayer at (p, z) = (0, — 4.5) (left curve, 
blue), (0,-2.5) (middle curve, magenta), and (0,-0.5) 
nm (right curve, red). The lipid bilayer extends from 
z = — 5toz = nm and is bounded by thin (black) 
vertical lines in Figs. [1] & [2j The cytosol lies below z = 
— 5 nm. The relative permittivities were taken to be 
Eu) = Ec = 80 and ei = 2. 

Figs. [1] & [5] show that the potentials fall off sharply 
as they cross the lipid bilayer. The reason for this ef- 
fect is a build-up of bound charge in the water near the 
lipid bilayer due to the high electric permittivities of the 
extracellular environment e^, and of the cytosol ec- 

In sections I Villi & IIXI I will extend this derivation 
to the case of several dielectric slabs. This generaliza- 
tion will allow us to add two layers of high-permittivity 
dielectric that will represent the head groups of the phos- 
pholipids. 



III. A SURFACE CHARGE ON A MEMBRANE 

As a pedagogical application of the formulas of the 
preceding section, let us consider a uniform charge den- 
sity a on the surface of a lipid bilayer of thickness t. To 
avoid minus signs, I will put the charge density on the 
extracellular leaflet. After doing the computation, I will 
translate the result to the case of phosphatidylserinc on 
the cytosolic leaflet. 

Maxwell's jump equation ([2]) tells us that the electric 
displacement Dm in the water differs by a from its value 
Di in the lipid, which in turn is the same as its value Dc 
in the cytosol. So we have two equations = D( + a 
and Di — Dc for three unknowns. We can use the elec- 
trostatic potentials of section|lT]to resolve this ambiguity. 

The electrostatic potential V^{p,z) due to a point 
charge at (0, 0, h) as given by ([9|) is equal to the electro- 
static potential V^{0, z) due to a point charge at (p, h). 
Thus setting the height /i in © equal to zero, and differ- 
entiating with respect to z, we find for the z-componcnt 
of the electric field at (0, 0, z) due to a charge q at (p, 0) 



h=0 



E,iO,z) 



47re^ 



{l+p)z 



(24) 



(25) 



{pp')"~^{z + 2nt) 

[p2 + {z + 27li)2]3/2 
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in which r = y7?+z^ and z > 0. Replacing the charge 
qhy a 2'Kpdp and integrating over p from p = to p = c», 
we have 



27rCT 
47re,„ 



1 + p -p^) ^{pp \ 

n=0 



(26) 



The dependence upon the variables z and t has dropped 
out. Doing the sum and using the definitions ((61) of p and 
p' , we get 



a (1+p)(1-pO ^ _ 
2e^, 1 - pp' Cu, 



PP' 



(27) 



The field points in the z direction. 

In the limiting case in which the electric permittivity 
of the extracellular medium is the same as that of the 
cytosol, e.u, = Ec, the electric field is i?w(cr) = ajle^j. 

We may apply similar reasoning to the formula \\\\ 
for the potential V'^{p^ z) in the cytosol due to a charge 
q in the extracellular water at a height h above the mem- 
brane. If we keep in mind that the quantity z — 2nt — h is 
negative, then we find for the electric field in the cytosol 
of a surface-charge density tr at /i = 



(28) 



Since there is no surface charge between the cytosol 
and the membrane, it follows from Maxwell's jump equa- 
tion @ that ~ €£ Eg = Dc ~ Cc Ec, and so that the 
electric field in the membrane is proportional to that in 
the cytosol 



Our formula 



E,{a) = -E,{a). 

for Ec{(t) now gives Ei{a) as 
ere,. 



Eeia) 



(29) 



(30) 



The jump in the displacement D across the layer of 
surface charge is 



D,„ -Df = 



[ew + Ec) 



(31) 



in agreement with Maxwell's equation ([2]). 

If the layer of surface charge of density, like that of 
phosphatidylserine, lies on the cytosolic leafiet at z = — t, 
then the electric fields are 



Eeia) 
Ecia) 



U{^m + ec) 
a 



(32) 



IV. IMAGE CHARGES 

As noted by Wagner , a charge q in or near a mem- 
brane polarizes the membrane and the surrounding wa- 
ter. The potential formulas l|9HTT|). (ITsHrfj), and (|2TH23l) 
represent these bound charges as infinitely many mirror 
charges. The mirror charges affect the behavior of ions 
near an interface between two dielectrics in ways that 
mean-field theories can't describe. 

For instance, a charge q in the lipid bilayer induces 
mirror charges in the cytosol and in the extracellular en- 
vironment. These induced charges are of opposite sign, 
and they attract the charge q in the lipid membrane. We 
can be more precise about this attraction if in the for- 
mula (fT6|) for V^^(p, z). we use V/(z) to represent the 
self-potential V^i^^z) without the z-independent, infi- 
nite, n = term of the first sum 



ln(l — pp') 



(33) 



P jpp')" 
^ \2z~2nt\ 



E 



p' wr 



\2z + 2{n+l)t\ 



Keeping only the first term in each sum and using C for 
the constant log term, we recognize two image charges 



q 



C 



P 
\2z\ 



P' 



\2z + 2t\ 



(34) 



familiar from freshman physics. They attract the charge 
q no matter what its sign. Water is better than lipid at 
attracting charges. 

Similarly, a charge q in the extracellular water induces 
a mirror charge in the lipid and others in the cytosol. 
The mirror charge in the lipid is of the same sign and, 
being closer, repels the charge q. We can describe this 
repulsion in terms of the formula for V^{p, z) if we 
use V^{z) to mean V^{0,z) without the z-independent, 
infinite term l/r 



47re„ 



P 
\2z\ 



pn-lp'n 



^'^i ^Ti |2^ + 2«t| 



E 



(35) 



The first term is the potential of the textbook mirror 
charge 



p 



47re,„ 12 z| 



(36) 



An ion of charge q in this potential has an energy pro- 
portional to q^ p, which is positive for both cations and 
anions. A lipid membrane therefore repels both cations 
and anions; the water attracts the ion more than the lipid 
does. In a mean-field theory, such as unpatched Poisson- 
Boltzmann theory, every particle responds to the same 
potential V{x), so the force qE{x) = —q'W{x) is pro- 
portional to the charge q of the ion and therefore must 
be opposite for cations and anions. Mean-field theories 
can't describe why a lipid membrane repels both cations 
and anions. 
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Ion Distributions near Charged Membrane 
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FIG. 3. (Color online) Monte Carlo predictions for the rel- 
ative concentrations of potassium pk{z)/Pk (solid, red) and 
chloride pci{z) / pci (dashed, blue) ions at a distance z (nm) 
from the charged cytosolic leaflet of a lipid bilayer are plot- 
ted along with the Gouy-Chapman predictions H40|) for the 
normalized potassium Kgc{z) (dot-dash,red) and chloride 
CIgc{z) (dots, blue) concentrations. 



V. IONS NEAR A CHARGED MEMBRANE 



FIG. 4. (Color online) Monte Carlo predictions for the av- 
erage total electrostatic potential V{z) (V) at a distance z 
(nm) from the charged cytosolic leaflet of a lipid bilayer as 
felt by a (solid, red) and by a CI" (dashed, blue) ion. 
The potential is that due to the phosphatidylserines of the 
cytosolic leaflet, the ions of the cytosol, and the polarization 
induced by the K'^ ion or by the CI" ion. The dot-dash ma- 
genta curve is the average potential felt by the ion without 
the self-potential H35p that represents the polarization the ion 
induces. 



One can use the formulas of Sec. |TT] in a Monte Carlo 
code to compute the distribution of salt ions near a 
charged membrane. Here I present the result of such 
a simulation of the distributions of potassium and chlo- 
ride ions near a membrane of a vesicle whose inner leaflet 
contains phosphatidylserine (PS) at a level of 4 percent, 
vifhich is about that of the plasma membrane of a liver 
cell. 

In the simulation, I let the potassium and chloride ions 
move according to a Metropolis algorithm within a box 
whose width and length were 50 nm and whose height 
was 10 nm. I took the potassium concentration to be 150 
mM so as to allow for a 10 mM concentration of sodium 
ions. The box contained 2258 K'^ ions. The bottom of 
the box was covered by a uniform negative surface charge 
density whose total charge was — 143 |e| corresponding 
to 143 phosphatidylserines at a molar density of 4 %. I 
used 2115 Cl~ ions to make the whole system neutral; 
these chloride ions played the role of the whole ensemble 
of anionic cell constituents. 

To mitigate edge effects, I surrounded the box with 
eight identical boxes into which I mirrored all 4373 ions. 
So there were 39,357 ions in nine identical boxes. All 
the boxes had the same uniform surface-charge density 
due to the presence of the PSs at a level of 4%. To 
strictly enforce periodic boundary conditions, one should 



use Ewald sums [10|, but I did not do this because they 
would have slowed the code down and because liquids 
are not crystals. Since the maximum step size in the z- 
direction was only I A; the error due to using only eight 
boxes of mirrored ions altered the energy difference Ai? 
of a Monte Carlo move by less than 0.064 fcT, usually 
about 0.016 /cT, far less than the thermal noise. 

In my Monte Carlo code [llj, I used the constant 
(P7)) for the electric field Ez (cr) of the surface charge of 
phosphatidylserines, the series ((35)) for the self-potential 
V^{z) arising from the response of the bound charge to 
an ion of charge g, and the sum © for the electrostatic 
potential V^{p, z) due to a charge in the extracellular wa- 
ter near a membrane. In this way, I took exact account 
of the fields of the charges of the problem while treating 
the neutral molecules of the extracellular environment, 
the membrane, and the cytosol as bulk dielectric media. 
The simulations consisted of eight separate runs in which 
23,000 sweeps were allowed for thermalization. Four of 
the runs collected data for an additional 50,000 sweeps; 
the other four for an additional 9,000 sweeps. 

In Poisson-Boltzmann theory p^ . all charges respond 
to a common potential V that obeys Poisson's equation 
with a charge density that respects the Boltzmann dis- 
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Relevance of Debye Layer for a Ion 
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FIG. 5. (Color online) For a potassium ion, the sum V^{z) + 
V„{z) of the electrostatic potential V^{z) due to the induced 
mirror charges (|35|) and that Vcr{z) = — a \z\/ {e^ + fc) due to 
the electric field (|32[) of a lipid bilayer whose cytosolic leaflet is 
charged to a phosphatidylserine mole percent of 4% (red, solid 
curve) or 20% (blue, dot-long-dash curve) is plotted against 
the distance \z\ (nm) from the leaflet. The uncorrected linear 
potentials Va{z) of the two surface-charge densities, 4% (blue 
dots) and 20% (red dashes), appear as straight lines. 



tribution 



lAV 



Ppf + 2^ Pfo^ e 



-QiV/kT 



(37) 



Here V is the electrostatic potential, ppf is a prescribed 
distribution of free charge, qi is the charge of species i 
and pfoi is its free-charge density where V vanishes p^ . 
This non-linear equation is hard to solve except in one- 
dimensional problems where the Gouy- Chapman solu- 
tion [11 is available [l2|. In the present context, that 
solutioir for the potential V{z) is 



where V =^ 0) arc then 

i^ccW = e-^^W/'=^ and CIgc{z) = e'^^'^^y'^^ . 

(40) 

In Fig.|3l I have plotted my Monte Carlo predictions for 
the relative concentrations of potassium iorrs pk{z)/pk 
(red solid curve) and of chloride ions pci{z)/pci (blue 
dashed curve) as functions of the distance z from the 
charged membrane. The Gouy-Chapman predictions 
(Uni) for the normalized potassium Kgc{z) (red, dot- 
dash) and chloride CIgc{z) (blue, dots) concentrations 
also are plotted there. The solid if"*" and dashed Cl~ 
Monte Carlo concentrations correctly drop sharply for 
respectively z < 1 and z < 2 nm due to the repulsion by 
the induced image charges as discussed in Section IIVI 
These concentrations are much lower than the Gouy- 
Chapman predictions. The Gouy-Chapman-Poisson- 
Boltzmann potassium concentration actually rises mono- 
tonically when the is less than 4 nm from the mem- 
brane. (The behavior of both pk{z)/pk and pci{z)/pci 
for z > 9 nm is an artifact due to the absence of ions at 
z > 10 nm in the simulation.) 

In Fig. m I have plotted my Monte Carlo predictions 
for the average total potential that a if + ion (solid, red) 
or a Cl~ ion (dashed, blue) feels due to the surface charge 
of the phosphatidylserines, to the other ions, and to its 
polarization of the three dielectrics. The induced bound 
charge sharply raises the potential felt by the potassium 
ion for 2 < 1 nm and lowers that felt by a chloride ion 
for z < 2 nm. The dot-dash magenta curve represents 
the potential due to the surface-charge layer of PSs and 
to all the ions (including the bound charges they induce) 
but without the image charges of equation ([55]) . These 
three potentials differ significantly over the whole range 
in which the Gouy-Chapman concentrations differ from 
their bulk values. (The dip in the three potentials for 
z > 9 nm is an artifact due to the absence of ions at 
z > 10 nm.) 

In the related papers 
potential. 



d, [ill, I neglected the sclf- 



VI. VALIDITY OF THE DEBYE LAYER 



V{z) = - — In 



\ — g-(2+zo)/A 



(38) 



in which A = 1/ s/SniBCoo, the Bjerrum length is is ~ 
/AireiukT, and Coo is the bulk ion concentration (taken 
to be the same for potassium and chloride). If a is half 
the absolute value of the surface-charge density of the 
phosphatidylserines, then the offset is 



zo = A In 



27r£BA 



- (l + ^1 + (2^^BAa)2) 



(39) 



The Gouy-Chapman formulas for the concentrations of 
the potassium and chloride ions (normalized to unity 



We have seen in the last two sections that mean-field 
theory cannot account for the behavior of ions near an 
interface between two dielectrics with two very different 
permittivities. Does this mean that the usual interpreta- 
tion of the Debye layer is incorrect? 

The answer depends upon the difference between 
the two permittivities and upon the magnitude of the 
surface-charge density. An image charge is proportioiral 
to the ratio of the difference of the two permittivities 
to their sum. So if the dielectrics have similar permit- 
tivities, then the induced charges will be weak, and the 
image-charge correction to a mean-field Debye layer will 
be small. Similarly, a high surface-charge density will 
dominate the field due to the induced bound charges un- 
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less the ions are very close to the interface. Thus the 
validity of the Debye layer depends on the relative mag- 
nitudes of the self-potential ~ qp/S-Keyj\z\ due 
to the image-charge correction ([55]) and the potential 
Va{z) = —a \z\/[ew + fc) due to the electric field (|32)) 
of the phosphatidylserines. In Fig. [5l I have plotted for 
a ion both their sum V^{z) -\- V^[z) and the po- 
tential Va{z). For the surface-charge density cr of a 4 
mole-percent concentration of phosphatidylscrine, as on 
the cytosohc leaflet of a liver cell, the sum V^{z) + Va{z) 
(red solid curve) differs from the surface-charge potential 
Va{z) (red dashed straight line) for z < 2 nm. But for 
a higher mole percent of 20%, the relative difference be- 
tween the sum V^{z) + Vcr{z) (blue dot-long-dash curve) 
and Va{z) (blue dotted straight line) is somewhat less. 

The importance of the image-charge correction (|35p 
rises as the surface-charge density falls. It is therefore 
particularly important in the case of an uncharged mem- 
brane, such as the outer leaflet of the plasma membrane. 



VII. A ZWITTERION 

Let us consider a simple model of a zwitterionic 
molecule in salty water above an uncharged lipid bilayer. 
The toy zwitterion is just a point charge q' at r' and an- 
other q" at r". The charges arc separated by s = r' — r" 
which makes an angle 9 with the vertical z so that 
z- s ^ s cos 6 where s is the distance between the charges 



s = Vs 



The square of the horizontal distance between 
the charges is = {x' — x")'^ + {y' — y")^ = sin^ 9. 
The midpoint of the molecule is r = (r' + r")/2 and its 
mean height is z ~ z ■ r. The heights of the charges are 
z' = z + is cos and z" ~ z ~ ^scos9. 

The electrostatic energy E'{z') of the interaction of 
the point charge q' with the polarization it induces is 



E'{z') = q''^u{z') = q'^u{z + cos 9) in which u{z) is the 
infinite sum (j35|) of image charges 



u(z) 



47re„ 



P 
\2z\ 



T — 



n=l 



2z + 2nt\ 



(41) 



Similarly, the function u(z) gives the energy of the inter- 
action of the other point charge q" with the polarization 
it induces as E"{z") = q"'^u{z") = q"^u{z - ^scos9). 

To find the energy E{r',r") of the charge q' in the 
fuU potential of the charge q" and the energy E{r",r') 
of the charge q" in the full potential of the charge 
q', we use our formula ([9]) for the potential V^{p,z). 
To compute E{r',r"), we imagine the charge q" to be 
at {x" ,y" , z") = (0,0, /i) and the charge q' to be at 
(p, 2,0) in cylindrical coordinates. Then E{r',r") is 
q'q"V^{p,z)/q where 



P 



q 4:new is ^/J^HVT^ 



(42) 



p'{i-p')j: 



{pp 



l\n-l 



\ + (z' + 2nt + z"Y 



> 
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FIG. 6. (Color online) The energy (|44l ) of a zwitterion in 
salty water a distance z from a lipid slab. The four lower, 
solid curves describe a molecule consisting of a point charge 
q = |e| separated by 0.5 nm from a charge q" = — |e|; the four 
upper, dashed curves are for a molecule consisting of a point 
charge q' = 2|e| separated by 2 nm from a charge q" — — |e|. 
Within each quartet the molecules from right to left make 
angles of 0, 30, 60, and 90 degrees with the vertical. 



This interaction is unchanged when we interchange the 
locations of the two charges, and z' + z" ~ 2z. So 
E{r',r") = E{r",r') = q'q"v{z,9,s) where 



v{z,9,s) 




P 



-p'{i-p')y: 



s2 sin^ 9 + (2z)2 
(pp')"-i 



(43) 



n=i ^s2 sin^ 9 + {2z + 2nt)^ 
The total electrostatic energy of the molecule is then 

E{z, 9) = q''^u{z + is cos 9) + q"'^u{z - ^s cos ( 



2 

+ 2q' q"v{z,9,s). 



(44) 



In Fig. [51 I plot the energy (jU ) of the model zwit- 
terion when it is z nm away from a plasma membrane 
considered as a lipid slab with salty water on both sides. 
The four lower, solid curves arc for a point charge q' = |e| 
separated by 0.5 nm from a charge q" = — |e|; the four 
upper, dashed curves are for a point charge q' = 2|e| sep- 
arated by 2 nm from a charge q" ~ — |e|. The lipid slab 
repels the molecules. 

One may find the energy of a more realistic zwitterion 
by integrating the two-charge formula (j44p over a suitable 
charge distribution. 
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VIII. SEVERAL DIELECTRIC LAYERS 



then 



We may extend our derivation of the electric poten- 
tial of a charge in a material of three dielectrics to the 
case of several dielectrics. One biological application is 
to a phospholipid bilayer considered as a layer of lipids 
bounded by two thin layers of head groups. These three 
layers with the extracellular water and the cytosol pose 
a five-dielectric problem. 

Let us first consider the case of four dielectric layers 
with the charge q in the first layer of permittivity ew in 
the region z > 0. Instead of the three potentials (jlj, we 
have four 



dk Jo{kp) 



-k\z-h\ 



u{k) e 



-kz 



dk Jo{kp) [mi(fc)e'' 
dkJa{kp) [m2(fc)e 



dk.h{kp)d{k) e^' 



kz 



(45) 



in which the first internal layer of permittivity ti fills the 
region —ti < z < Q while the second internal layer of 
permittivity £2 fills the region — <i — t2 < z < — ti. The 
constraints ([T] & [2]) give six equations 

"M + /l - u = /5 
eiTOi - ei/i + e^u = e.„/3 
TOi + 2/1/1 - "i2 - yi/2 = 
eiyi/i - <^2-m2 + e2yi/2 = 



eimi 



m2 + 2y2./2 - d = 
e2TO2 - e22/2./2 - Ecrf = (46) 

in which yi{k) = exp(2fcfi) and y2{k) = exp(2fc(<i -l-t2)), 
while as in ([S]) the parameter /3(fc) ~ q exp{—kh) / Anein 
represents the charge at z = ft, > 0. The functions mi(fc) 
and /i(fc) determine 7712 (fc) and f2{k) as 



m2 = [(£2 + ei)mi + (£2 - ei)?/i/i]/2e2 
/2 = [(£2 - £i)mi/yi + (£2 + ei)/i]/2e2. 



(47) 



We now address the problem of a charge q in a semi- 
infinite region of permittivity e^, at a height h above n 
internal layers of permittivity et and thickness U which 
in turn are above a semi-infinite region of permittivity 
Ec- In the ith internal layer, the potential is 



vrip,z) 



dkJoikp) m,(fc)e*^ + /,(fc)e 



-kzl 



(48) 



while the first and fourth equations of the set p5|) de- 
scribe the potentials in the two semi-infinite regions. The 
functions mi{k) and fi{k) determine those mi+i{k) and 
/,;+i(fc) of the next layer by matrix multiplication. If 



Pi 



and 



Ci+l 



Cj+l 



(49) 



fi+i 



1 



Q+i + Ci (e,;+i - ei)yi\ f m, 



2e«+i 



(50) 



in which yo = I and yi = exp[2fc(ii -I- • • • -I- U)] for i > 0. 
Let us set mo{k) = f3{k) and fo{k) = u{k) as well as 
TO„+i(fc) — d{k) and /„+i(fc) = 0. In these formulas, 
^n+i is 6c, and 69 is 6^^,, not the permittivity of the vac- 
uum. If is the product of i' + 1 of the matrices (|50p 



then we have 



i=0 



fi+i 



PiVi 



,+i \Pi/yi 1 



(51) 



(52) 



Setting i ^ n gives us 

u = -pEi-^/Ei-^ 



^22 



(53) 



IX. THREE DIELECTRIC LAYERS 

A phospholipid bilayer consists of a layer of phosphate 
head groups, a (double) layer of lipids, and a second layer 
of phosphate head groups. In this section, we will apply 
the formulas of section IVIIII to the problem of a charge q 
at a height h above such a bilayer. There are now three 
slabs and two semi-infinite regions. We must compute 
£:(3). I get 



E. 



(3) 

21 



Po 



El 

yi 



El 

y2 



PQPiP2yi 
2/2 



P3 I 1 1 

— 1 +popiyi 
2/3 



' P0P2y2 



piP2y2 \ 
yi ) 



(54) 



and 



E. 



(3) 
22 



POPI , P0P2 , PlP2yi 



yi 

P3 ( 

— [pa- 



y2 



y2 



piyi +P22/2 



P0PlP2y2 \ 

yi J 



(55) 



Stern and Feller [14|, Nymeyer and Zhou [15|, and 
Baker [3l have estimated the relative electric permit- 
tivity of phospholipid membranes as being 1 from to 
10 A from the center, 4 from 10 to 15 A, 180 from 15 to 
20 A, 210 from 20 to 25 A, and like bulk water beyond 25 
A. I will approximate their results by using 2 from to 
15 A and 195 from 15 to 25 A and will take the relative 
permittivities of the cytosol and of the extra-cellular envi- 
ronment to be 80. These approximations greatly simplify 
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FIG. 7. (Color online) The potential V^[p,z) (V,[59| at the 
point (p, z) due to a charge g = |e| on the z-axis at (0, h) above 
a three-slab phospholipid bilayer that is neutral (lowest, solid, 
blue) or has a 4 mole-percent layer of phosphatidylserine on 
its cytosolic leaflet (upper, solid, red). Both p and h are 1 nm. 
The potential @ of a charge above a single neutral lipid slab 
without head groups is plotted as a dashed magenta curve to 
illustrate the effect of the layers of head groups. 



FIG. 8. (Color online) The potential V^{z) (V, EH) feh by 
a unit positive charge at a height z nm above a phospholipid 
membrane modeled as a lipid slab bounded by two thin polar 
slabs. The lowest (solid, blue) curve is for a neutral mem- 
brane; the upper (solid, red) curve is for a membrane whose 
cytosolic leaflet is negatively charged by phosphatidylserine at 
4 mole percent. Without the outer polar slabs, the potential 
psp is purely repulsive (dashed magenta). 



The key function u{k) then is by ([5^ the ratio 



our formulas JSH&ISn]) and imply that po = 0.418 = — ps 
and pi = —0.98 = —p2- The thicknesses of the layers are 



t 



3i, and is 

AOkt 



t, and so yi = e 



2kt 



t = 1 nm, t2 

y2 = e^*"'*; and ys — e^*^*^*. With these simplifications, 
our formulas ([M] & (55]) reduce to 



e: 



(3) 
21 



poe-^'^^l+popiie'"* 



2 ~6kt 

PoPie 



„8fet 



= Po+PoPie 
-Poe-''''- 



-4kt 



+ Pi(l+Po)(e 



-2kt 



-8kt\ 



(56) 



2 -6kt 

PoPie 



and 



El^f = l+popi{e 



~2kt 



'8kt\ 



-Poe-''"' {po+Pi{e''' - e''')-pople'''') 



(57) 



u{k) 



1 _^-kn 



Po +Pi(l +Po)(e 



2\(^-2kt „-ikt\ 



1 + 2popi{e 



-2kt 



-8kt 



)+plple' 



-Akt 



p2e-6fct_p2g-10/ct 



(58) 



We can use it and our formula (|45|) to write the poten- 
tial V^{p, z) at the point (p, z) due to a charge q on the 
z-axis at (0, h) above a three-slab phospholipid bilayer as 



dk Jo{kp) 



q 



-k\z-h\ 



uik) e 



-kz 



(59) 

Figure [7] plots this potential V^(p, z) for p = 1 nm and 
< z < 5 nm for the case of a unit positive charge 9 = |e| 
on the z-axis at height h — 1 nm above a three-slab 
phospholipid bilayer that is neutral (lowest curve, solid, 
blue) or has a 4 mole-percent layer of phosphatidylserine 
on its cytosolic leaflet (upper curve, solid, red). The 
electric field of the PSs is taken from to be [a) = 
+ Cc) in which a is negative. To illustrate the 
effect of the layers of head groups with very high electric 
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permittivity, I have replotted the potential ([9]) of Fig. [T] 
due to the same charge but above a single naked, neutral 
lipid slab as a dashed magenta curve. The head-group 
dipoles lower the potential V^{p, z) in their vicinity. 

We also can use the ratio u{k) to compute the self- 
interaction of a charge q with our three slab model of the 
phospholipid bilayer. The potential felt by the charge in 
the salty water due to all the image charges its presence 
induces in the three slabs and in the cytosol is 



would repel. If the membrane has phosphatidylserines on 
its cytosolic leaflet, then it strongly attracts positive ions 
that it otherwise would repel. Apparently the phosphate 
head groups facilitate many physiological processes, such 
as the docking of ligands and the translocation and en- 
docytosis of positive ions and cell-penetrating peptides. 



X. SUMMARY 



u{k) e-"^ dk 



(60) 



in which u{k) is the ratio ((58)) but with the height h of 
the charge q replaced by z. The potential felt by the 
charge q at z then is 



-2kz/t 



, 2 -4fc 

Po+PoPie 



Poe 



-lOk 



PoPie 



1 



-pMe-'" 



pIc'" 



lOfc 



(61) 



-f- 2popi[e 



-2k 



1 



p\e~ 



■6k 



dk 



in which the parameter t is one nm. 

Figure [S] plots the electric potential felt by an ion of 
charge |e| near a phospholipid bilayer modeled as a lipid 
slab bounded by two thin polar slabs. The lowest (solid, 
blue) curve is for a neutral membrane; the upper (solid, 
red) curve is for a membrane whose cytosolic leaflet is 
negatively charged by phosphatidylserine at 4 mole per- 
cent represented as in Fig.[71 The two thin layers of phos- 
phate head groups make the potential felt by an ion near 
a neutral phospholipid bilayer attractive rather than re- 
pulsive. Without the outer polar slabs, the potential ([35)) 
is purely repulsive (dashed magenta). The head groups 
attract to the membrane ions that a naked lipid slab 



I derived the electrostatic potential of a charge in or 
near a lipid bilayer in section HIl and used it in section Hill 
to compute the electric field of a uniformly charged mem- 
brane and in section IIVI to describe the effects of image 
charges. In section |Vl I used the results of sections IIIllIVI 
in a Monte Carlo computation of the distribution of ions 
near a charged membrane. I discussed the validity of the 
Debye layer in section I VI I and computed the energy of a 
zwitterion near a lipid slab in section lVlIl In section lVIII[ 
I calculated the potential of a charge near a membrane 
modeled as several dielectric layers of different permit- 
tivities between two different semi-infinite dielectrics. I 
used this analysis in section llXl to model a phospholipid 
bilayer as a lipid layer bounded by two layers of head 
groups of high electric permittivity. The phosphate head 
groups cause a neutral membrane to attract rather than 
to repel ions. 



ACKNOWLEDGMENTS 

I am grateful to Nathan Baker, Leonid Chernomordik, 
Charles Cherqui, David Dunlap, Scott Feller, Kamran 
Melikov, Michael Wilson, Adrian Parsegian, Sudhakar 
Prasad, Harry Stern, and David Waxman for helpful 
conversations; to Susan Atlas, Bernard Becker, Vaibhav 
Madhok, Samantha Schwartz, James Thomas, and Toby 
ToUey for useful comments; and to the two referees for 
constructive criticism. 



[1] E. M. Severs, P. Comfurius, D. W. Dekkers, and 
R. F. Zwaal, Biochim Biophys Acta 1439(3), 317 (1999); 
B. Alberts, A. Johnson, J. Lewis, M. Raff, K. Roberts, 
and P. Walter, "Molecular biology of the cell," (Garland 
Science, New York, 2002) pp. 587-593, 4th ed. 

[2] G. L. Gouy, J. Phys. 9, 457 (1910); D. L. Chapman, 
Philos. Mag. Ser. 6, 25(148), 475 (1913); C. Wagner, 
Phys. Zeits. 25, 474 (1924). 

[3] L. Onsager and N. N. T. Samaras, J. Chem. Phys. 2, 528 
(1934). 

[4] A. Parsegian, Nature 221, 844 (1969); D. Walz, E. Bam- 
berg, and P. Lauger, Biophys. J. 9(9), 1150 (1969); 
B. Neumcke and P. Lauger, ibid. 9(9), 1160 (1969); 
B. Neumcke, D. Walz, and P. Lauger, ibid. 10(2), 172 
(1970); R. H. Brown, Prog. Biophys. Mol. Biol. 28, 341 
(1974); A. Parsegian, Ann. NY Acad. Sci. 264, 161 



(1975); B. Jonsson, H. Wennerstrom, and B. Halle, J. 
Phys. Chem. 84, 2179 (1980); S. McLaughlin, Ann. Rev. 
Biophys. Biophys. Chem. 18, 113 (1989); R. R. Netz, 
Eur. Phys. J. e'3, 131 (2000); Y. Levin and J. E. Flores- 
Mena, Europhys. Lett. 56, 187 (2001); R. Allen and J.-P. 
Hansen, Mol. Phys. 101(11), 1575 (2003); Y. Levin, Eu- 
rophys. Lett. 76(1), 163 (2006); Y. S. Jho, G. G. Park, 
C. S. Chang, P. A. Pincus, and M. W. Kim, Phys. Rev. 
E 76, 011920 (2007); A. P. Shreve, M. C. Rowland, A. R. 
Sapuri-Butti, T. W. Allen, and A. N. Parikh, Langmuir 
24, 13250 (2008); D. Frydel, J. Chem. Phys. 134, 234704 
(2011); N. Shimokawa, S. Komura, and D. Andelman, 
Phys. Rev. E 84 (2011). 
[5] A. Bakhshandeh, A. P. dos Santos, and Y. Levin, 
jPhys. Rev. Lett. 1077107801 (2011)| 



12 



[6] G. V. Miloshevsky, A. Hassanein, M. B. Partenskii, and 
P. C. Jordan, J. Chem. Phys. 132(23), 234707 (2010). 

[7] M. A. Wilson and A. PohoriUe, JAGS 118(28), 6580 
(1996); S. Dorairaj and T. W. Allen, PNAS 104(12), 
4943 (2007). 

[8] J. Schwinger, L. Deraad, K. A. Milton, and W.-y. Tsai, 
"Classical Electrodynamics," (Westview Press, 1998) 

Ghap. 16. 

[9] K. 



Phys. Biol. 7, 016001(14pp) (2010) 



E. GahiU, 
arXiv:0810.2358. 
[10] P. P. Ewald, Ann. Phys. 369, 253 (1921). 



[11] K. E. GahiU, "Fortran-90 codes for the monte 
carlo simulations of this paper," (2011), 
http://bio.phys.unm.edu/membElec/index.html. 

[12] P. Nelson, Biological Physics, updated first ed. (W. H. 
Freeman, New York Gity, 2008). 

[13] K. E. GahiU, |IET Syst. Biol. 4(6), 367 (2010) 
arXiv: 1009.3606. 

[14] H. A. Stern and S. E. Feller, J. Ghem. Phys. 118(7), 
3401 (2003). 

[15] H. Nymeyer and H.-X. Zhou, Biophys. J. 94, 1185 (2008). 
[16] N. Baker, (2011), e-maU. 



